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It is proposed in this paper that the total lunar echo at any instant of time is comprised 
of two components, which are the result of scattering from two general types of terrain fea- 
tures. On the basis of this model, the average theoretical backscattered power flux at the 
receiver is derived as a function of time for the case of short pulse transmission. The radar 
cross section of the moon is also calculated. Several comparisons are made between the 
theory and experimental results obtained with a variety of equipment parameters. 



1. Introduction 

More than a quarter century ago experimenters 
initiated efforts to establish radar contact with the 
moon. It was not until 1946 that radar equipment 
had been improved to the point where positive re- 
sults were obtainable. Since that time the problem 
of radar reflections from the moon has been subject 
to extensive investigation by several experimental 
groups and theoretical workers botli here and 
abroad. 

Among the hrst theoreticians to report their find- 
ings in this general area were Grieg, Metzger, and 
Waer [1948]. Although only continuous wave trans- 
mission was considered, the total transmission loss 
over the earth-moon-earth path was derived for both 
the smooth and the perfectly diffuse lunar surface 
models. 

During the past few years several experhnental 
groups have looked for moon echoes using pulse 
transmission at a variety of frequencies. Among 
these investigations are those reported by Trexler 
[1958], Yaplee et al. [1958], Hughes [1960], Leada- 
brand et al. [1960], and Pettengill [1960, 1961]. 
Adaptations of neither the smooth nor the diffuse 
moon models could account for the shape of the 
returns using short transmitted pulses. This dis- 
parity is clearly illustrated in figure 1 which presents 
experimental results obtained by Pettengill [1961] 
at a freciuency of 440 Mc/s using a transmitted pulse 
duration (r) of 100 jitsec. If the moon scattered 
radiation like a smooth sphere then the echo would 
consist almost entirely of a pulse return with a dura- 
tion of the order of 100 jusec. In the other extreme, if 
the lunar surface were a completely diffuse reflector 
at radar wavelengths then the initial rapid rise of the 
return would be followed by a much slower decay 
than is actually observed. Theoreticians were 
thereby stimulated to begin their studies of the limar 
reflection problem anew. 
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Senior and Siegel [1959] first proposed what might 
be termed a ^^phenomenological" theory in which it 
was suggested tliat the lunar echo consisted of a 
small number (4 or 5) of individual returns, each of 
which was associated with a single reflecting area. 
However, other workers (e.g., Hughes [1960]) felt 
that experimental data obtained with very short 
pulses failed to support this view. In response to 
such criticism, Senior and Siegel [1960] later increased 
their estimate of the niunber of scattering areas to 25 
or 30, all of which were located within a depth of 
about 50 km from the subterrestrial point of the 
moon. However, a revision of this nature fails to 
resolve two rather serious objections. First, the 
theory as it is formulated does not provide a func- 
tional formalism which pei'mits quantitative compar- 
ison with experimental results obtained using a 
variety of piflse lengths. Secondly, on the basis of 
our ciUTent knowledge of the moon's surface, it seems 
quite unlikely that a number of large relatively 
smooth ''spherical siu'faces" are concentrated only in 
the immediate vicinity of the subterrestrial point, 
while the remainder of the scattering regions over the 
lunar surface are ''craterlike. ^' In any event, the 
theory of Senior and Siegel lias not received wide 
acceptance and most recent work has used a statis- 
tical approach. 

The basic premise of the statistical theories is that 
at any instant of time the amplitude of the lunar echo 
represents the net efl'ect of scattering from a large 
number of comparatively extensive reflecting areas 
randomly oriented and distributed over the siuface. 
Several workers have reported theoretical results in 
this area, notably Hargreaves [1959], Daniels [1960a, 
1960b], Brown [I960], and Hagfors [1961]. These 
theories need not be reviewed here since a brief 
summary of each is available in a recent survey paper 
by Evans [1961]. 

Section 2 of this paper is devoted to a discussion of 
lunar scattering mechanisms in which it is proposed 
that the total echo at any instant of time is the result 
of reflection from two general types of terrain 
features. One of these is described by a statistical 
model which is discussed in some detail in section 3. 
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The radar cross section of the moon is calculated in 
section 4. Sections 5 and 6 are devoted to deriva- 
tions of the temporal distribution of lunar echo power 
in the case of short pulse transmission. Finally, 
several comparisons between theory and experiment 
are given in section 7. 

2. Discussion of Scattering Mechanisms 

The experimentally obtained echo amplitude dis- 
tribution appearing in figure 1 leads one to suspect 
that perhaps more than a single scattering mechanism 
is responsible for the observed decay. Clearly, the 
initial part of the return is suggestive of a specular 
type of scattering process rather than a diffuse type. 
On the other hand, the subsequent portion of the echo 
reveals just the opposite characteristic. Thus, it 
seems appropriate to regard the total signal at each 
instant of time as the combined eifect of two scatter- 
ing components of different types which, for con- 
venience, will be referred to hereafter as the '^quasi- 
specular" and the '^diffuse" components, respectively. 

The quasi-specular return differs from that which 
would be observed from a smooth moon in that there 
exists an extended period of signal decay at times 
greater than the original pulse length. As pointed 
out by Evans [1961], the results of experiments using 
different pulse lengths suggest that during this period 
of decay the echo amplitude represents the net effect 
of reflection from a large number of randomly 
distributed scattering areas. This observation leads 
quite naturally to the so-called ' 'statistical" model 
for quasi-specular scattering from the lunar surface. 
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Figure 1. Observed echo power as a function of time for a 
pulse lengtKjof 100 ixsec. 



According to this model, the moon is characterized 
by random irregularities in surface height h. Further- 
more, these irregularities are presumed to be large 
in the sense that their characteristic dimensions are 
somewhat greater than a few meters and that other- 
wise the surface is smooth, with the result that in- 
coming radiation is scattered locally in a specular 
fashion. From the mathematical point of view, it 
is supposed that h can be treated as a continuous 
random function of position on the lunar surface and 
that the statistics of the function h are known, apart 
from characteristic parameters which are to be 
determined from experiment. 

The experimentally obtained echo decay rate is 
determined as an average over a large number of 
returns. Since the rough lunar surface presents 
a slightl}^ different aspect from one pulse to the 
next as a result of lunar libration, one might think 
of each experimentally obtained return (or set of 
returns, depending upon the prf) as representing 
the echo obtained from a member of an ensemble 
of moons characterized by different surface height 
functions h. In an analogous fashion, the average 
amplitude distribution of the return predicted on 
the basis of the statistical model is determined by 
obtaining an ensemble average of the expression for 
returned power. The average echo distribution 
will resemble specular reflection much more than 
diffuse reflection so long as the deviation of the 
actual surface from the mean is not large and the 
average surface height-to-separation ratio character- 
istic of the irregularities is fairly small. On the 
other hand, it seems quite unlikely that such a 
scattering mechanism could explain the slow, low- 
level portion of the decay which follows the initial 
rapid amplitude decrease. 

Although the statistical model embodies those 
ph3^sical characteristics of the surface which play a 
dominant role in backscattering near normal inci- 
dence, it is incomplete in the sense that it does not 
include the existence of small irregular surface 
features, such as rocks and jagged facets, which 
will tend to scatter radiation more isotropically. 
These latter irregularities will contribute relatively 
little to the amplitude of the signal scattered from 
the vicinity of the subterrestrial point of the lunar 
surface. On the other hand, in the case of a short 
transmitted pulse the echo amplitude at times 
considerably greater than the pulse length is due 
to returns from portions of the moon's surface 
closer to the limb. In this region reflections from 
small, less directive scattering facets w^ill assume 
greater importance than quasi-specular reflection 
from the slow undulatory component of surface 
roughness. Therefore, it is natural to expect that 
scattering from these smaller irregular facets will 
account for the slow decay of echo amplitude which 
follows the rapid decrease from the peak of the 
return due to quasi-specular reflection. 

In view of the fact that the quasi-specular and 
diffuse reflection mechanisms operate essentially 
independently of each other, the total lunar echo 
power at each instant of time can be considered to 
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be tlir sum of the scattered powers associated with 
the individual components. Consequently, the 
power amplitude due to these different scattering 
mechanisms can be derived independently and the 
results can subsequently be added to give the com- 
posite lunar echo. Since the maximum response 
from diffuse scattering is 15 to 20 db below that due 
to quasi-specular reflection, it is apparent that ordy 
the latter component need be considered in estimat- 
ing the total return in the case of continuous wave 
transmission. 

3. Statistical Model of the Lunar Surface 

The center of the moon is taken to be coincident 
with the origin of a Cartesian coordinate system as 
illustrated in figure 2. The transmitter-receiver is 
located on the positive 3:3-axis at a distance R from 
the origin. 

The surface of the moon is treated as tliough it 
were randondy I'ough. Thus, if the mean surface 
is assumed to be a sphere of radius a, then the local 
radius r at any point X = (0^1,^-2,^3) on the true surface 
can be expressed as r = a-\-h, where the local surface 
height h is a random function of position. In ex- 
pressions involving the surface configuration it is 
often convenient to employ spherical coordinates, so 
that the locus of the point x can be represented by 
the pai'anu^l I'ic syst(Mn 



Xi = r sin 6 cos 
X2=r sin 6 sin 
Xs=r' cos 6 



(3-1) 



with r = a+h(d,4)). The local normal vector at x 
can be calculated from the expression 



n=[J,e,+J2e2+Jse,]/[Ji+Jl+JI]^ 






(3-2) 



where e7,:(z= 1,2,8) are the Jacobians of the system 
(3-1). The ti are unit vectors in the three mutually 
orthogonal directions. The local surface area ele- 
ment dA is given by 



dA=[Jl+Jl+JI]^ddd<l>, 



(3-3) 



It is consistent with the hypothesis of locally spec- 
ular scattering to assiune that those surface rough- 
nesses which determine the |)roj)erties of the ({uasi- 
specular component are characterized by dhnensions 
which are greater than several radar wavelengths. 
In the sequel we shall denote by a the standard 
deviation of the local surface height, and 8 will 
denote the great circle distance separating two 
points on the mean lunar surface beyond which local 
surface heights are virtually uncorrelated.^ In 
view of the fact that as a radar reflector the near 
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Figure 2. Geometry of the Iniuir scattering problem. 



surface of the moon ap])ears to bo much more smooth 
than rough, it is ant ici patent that the vahie of a/d 
will be at least an order of magnitude less than unity. 
This remark follows from the observation that a/d is, 
in a sense, a measure of the rms surface height-to- 
separation ratio. Effects such as shadowing and 
multiple reflections will therefore be of secondary 
importance near the subterrestrial point of the moon 
and need not be considered here. 

An analysis based on a statistical model of the 
present sort becomes prohibitively complicated 
when an attempt is made to take into account 
explicity the possible electromagnetic properties 
of the moon^s crust. On the other hand, the statisti- 
cal model proposed here implies a gently undulatory 
lunar surface, which is smooth over distances of a 
few radar wavelengths and which scatters radiation 
specularly in a local sense. This suggests that any 
given region of the moon's surface must contain 
areas which are favorably oriented to produce re- 
turns if it is to be at all eft'ectivc in contributing 
to the back scattered power. This leads to the 
assumption that the relative importance of any given 
region in contributing to the returned echo will be 
more strongly dependent upon the local geometry 
than upon the electromagnetic parameters of the 
surface material. In order to determine the in- 
fluence of local geometry upon the observed echo 
it has become standard practice in discussions of 
lunar scattering to assume that the far zone scattered 
fields can be calculated by taking the total flelds at 
the lunar surface to be the same as though the 
interior were a perfect dielectric with e^= oo . The 
scattered power level in the far zone is later adjusted 
by a multiplicative correction factor which is usually 
referred to as an effective power reflection coef- 
ficient. The value of this coeffici(*nt is determined 
by calculating the fractional power scattered by a 
semi-infinite dielectric slab with an (^ffcH'tive relative 
permittivity €r exposed to a normally incident plane 
wave. Since the surface areas providing the greatest 
contribution to the quasi-specular scattering com- 
ponent are those which are perpendicular to the 
direction of propagation of the incident field, this 
procedure is a good approximation. 

Since a dielectric with er = oo is indistinguishable 
from a perfect conductor the possibility of internal 
reflections is implicitly eliminated from considera- 
tion. This is consistent with reality since the lunar 
radius a is so much greater than the skin depth of 
rocky material at radar wavelengths. 
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The section which follows is devoted to the devel- 
opment of an approximate expression for the radar 
cross section of the moon on the basis of the scat- 
tering model just described. 

4. Radar Cross Section of the Moon 

At the conclusion of section 2 it was remarked 
that the total radar cross section of the moon may 
be calculated with a fair degree of accuracy by 
assuming that only the quasi-specular component 
contributes significantly to the backscattered power. 
This calculation will be performed using the statis- 
tical model for the lunar surface. 

In the case of continuous wave transmission we 
may represent the incident magnetic field as a plane 
wave propagating in the negative Xa-direction with a 
time dependence of the form exp {—io)t)\ 



lli^^iHi exp [—i{kx^+oit)Y 



(4-1) 



The scattered magnetic field can then be calculated 
from the integral of Stratton and Chu: 

X(nXH) + (n.H)v(7]rfA (4-2) 

In this expression, 6^ is the standard Green's function 
for the three dimensional Helmholtz equation and 
the electromagnetic vectors in the integrand are the 
total fields at the surface of the moon. The vector 
n is the outward local normal to the true surface at 
the point (0,0). Since reflections from the portion 
of the moon's surface in the shadow region are 
neglected, the surface A of integration may be taken 
to be the illuminated hemisphere. 

The several assumptions outlined in the previous 
section permit the use of the so-called ' 'tangent plane 
approximation." This, together with the idealized 
electric properties of the surface material, provides us 
with simple relations for the fields at the surface: 



and 



H = 2H,-2(Hrn)n 

nXE = 0. 



(4-3) 
(4-4) 



Following the substitution of (4-3) and (4-4) into 
the integrand of (4-2) and the specialization of the 
result to the far zone, the scattered magnetic field 
can then be written in the form 

tik "I 
Y-jn U^^exp [i{kR 

—oit)] j exp i—i2kxs)Jsddd(t). (4-5) 
The far zone scattered electric field can be written as 



Es=-toe2Hs, 



(4-6) 



where fo is the characteristic impedance of free space. 
The quantity of particular interest is the expected 
value of the Xs-component of the Poynting vector S 
at the position of the observer. In the case of fields 
varying harmonically in time, it is customary to sup- 
press the rapid fluctuations of S by calculating a local 
time-averaged value. This operation is certainly 
meaningful for continuous wave transmission when 
the surface of the scatter is characterized by irregu- 
larities considerably greater than a wavelength in 
dimension. Thus the quantity to be derived is the 
expected value of 



S=^ReE,XHr 



(4-7) 



Under the present set of assumptions the Poynting 
vector at the receiver is directed along the cc3-axis so 
that no confusion will arise if Sx^ is written simply as 

Equations (4-5) and (4-6) can be combined with 
(4-7) to give an expression for S in terms of the 
random surface height h. The expectation or 
stochastic average of the power flux is obtained 
through multiplication of S by an appropriate joint 
probability density function pihiM followed by 
integration over all hi and A2. This general method 
was first presented by Isakovich [1952] and later 
placed on a more rigorous basis by Hoffman [1955], 
both of whom treated reflections from a rough 
conducting plane. After an application of Fubini's 
Theorem this procedure in the present case leads to 
the relation ^ 

Re exp [—i2ka{cos ^1— cos ^2)] 

•(^{exp [—i2k(hi cos di—h2 cos ^2)] 

JuJz2]ddidd2d(i>4<t>2 . (4-8) 

In order to evaluate this last expression, a number of 
simplifications must be introduced. 

The quantity J3 is easily calculated from system 
(3-1), with the result 



'[ 



J->=r^ sin d I cos d-\ — ^- sin B 



'} 



(4-9) 



First we observe that h is generally expected to be 
an extremely small quantity compared with a, so 
that little error is incurred if r^ in the first factor is 
set equal to (j?. The second term inside the paren- 
theses is nothing more than the product of sin d and 
the local surface slope in the ^-direction. Current 
estimates of the rms value of lunar surface slopes 



3 Expectations are denoted by ^ { } to avoid confusion with the symbol £" which 
is reserved for the electric field. 
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place this fi^^iiro at about 10° or less. Lunar radar 
data ol)taine(l by using very short pulses (see fig. 1) 
indicate that n(\irly all of tlie power returned by the 
lunar hemisphere is scattered from a region in the 
vicinity of the subterrestrial point. It would appear 
that if the estimates of the lunar slopes are reasonably 
accurate, the second term in parentheses in (4-9) 
does not assume importance compared with the first 
except in those regions of the surface which con- 
tribute negligibly to the total radar return. This 
term will henceforth be omitted. 

As a result of the simplifications just outlined, Jz 
no longer involves the random variable A so that e/31 
and J32 may be removed from the braces. The ex- 
pectation which remains is simply the characteristic 
function appropriate to the joint probability density 
V^iM' If we assume this latter function to be the 
standard Gaussian bivariate distribution with vari- 
ance 0-^ and correlation function p, then (4-8) can be 
expressed in the form 

^^{S} =[§£??] [^ J Re [£ J^ exp[-i2^a(cos^i 

—cos (92)] • exp [— 2ArV(cos<9i 

— cos 6/2)^] exp [— 4/:V cos Qi cos ^2 (1 — p)] 

• cos (9i cos 0. sin B^ sin B-. dd^didftyM-i . (4-10) 

We first observe that the third exponential factor 
in the integrand of (4-10) is a rapidly varyhig func- 
tion when p is near unity. This suggests that we 
replace p by the first few terms of an appropriate 
series expansion. The correlation function p is near 
unity only when the points (^i,</)i) and (^2,</>2) on the 
mean surface are close to one another. Therefore, 
throughout that portion of the coordinates domain 
where the series expansion is no longer valid, the 
harmonic factor oscillates very rapidly in general, 
and tends to cancel out residual errors introduced by 
the expansion. 

The functional form of p is somewhat arbitrary, 
but in view of the fact that the mean surface is a 
sphere, it is natural to assume it to be a function of 
x/2, where X is the angle between the lines joining 
points (^i,0i) and {62,^2) with the center of the sphere. 
A large class of correlation functions, for example, is 
represented b}^ the relation 



p=p(sin x/2). 



(4-11) 



?^ If we exclude the possibiHty of a discontinuity in 
I slope at X=0, we can expand p in powers of sin x/2 

\ to obtain 

I 

j 

P==l-^|p''(0)|sin2x/2+... 

f ==l-i|p''(0)|[l-cos6'icos(92 

— sin ^1 sin (92 cos (0i— (/)2)] + . . ., (4-12) 



where the prime indicates differentiation with respect 
to sin x/2. 

It is now expedient to take further advantage of 
the fact that nearly all the power is scattered from 
the vicinity of the specular reflection point. If we 
put iCi^sin Oi in (4-10) then most of the radar signal 
is returned from the region of small Xi. Thus we 

shall set cos ^^=1—- x\ throughout the integrand 

and drop terms which are of an order higher than x\ 
where appropriate. Further simplification results 



is equal to the product of 



from the fact that ^p'^(O) 

{2a?l(T^) and the mean square value of the local sur- 
face slope and is much greater than unity except for 
completely negligible slopes. (It may be noted that 
in the case of a Gaussian correlation function this 
product is equal approximately to 4a^/5^) Thus we 
may neglect the terms Ba^x\ as compared with 



1 



Z:V|p''(0)|:rf,(t=l,2). With X^=P(r2|p''(0)| and v 
--2ka, the expectation of S takes the form 



Si 



'^'<l'^0[S.J''» 



L{k)}, (4-13) 



where 
L 



m2ir r2- 
I Jo 



exp |^i^(xi-^^) J exp |^-- (xj+j^) J 

• exp [X^XiXo cos (</)i— </)2)] dx{dxld<t)i(l<j)2. (4-14) 
If we now put }Xi= — -{iv — \^) and jlx2=^(^V+X^), 

then the result of integration over the azimuthal 
angles in (4-14) can be written as 



L{k) =47r^ exp [— (mi^i+M2^2)] h()<^XiX2)dx\dxli 

Jo Jo 

(4-15) 

where /o is the zero order modified Bessel function 
of the first kind. Since X^ is a number very much 
greater than unity, the Bessel function can be ex- 
panded asymptotically for Xi of the order of unity 
and greater to show that little error is incurred if 
the upper limits of integration in (4-15) are set equal 
to infinity. If /o is then replaced by its series ex- 
pansion and the summation is removed from the 
integral signs, then the integrations are easily per- 
formed to give as a result 



L{k)- 



4t^ 
M1M2 



m=0 L4M1M2J 



4^ 

''k^a^' 



(4-16) 
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Tims, in the present approximation the expected 
power flux at the receiver is independent of the pre- 
cise form of |p''(0)|. The expected value of S is 
given by 

^iS}-[^Hq[^J (4-17) 

When (4-17) is multiplied by an effective reflection 
coefficient a, standard radar formulas can be used 
to show that the effective radar cross section of 
the moon in the present approximation is 

as = a7ra\ (4-18) 

This is exactly the same result obtained by Grieg, 
Metzger, and Waer [1948] for a smooth moon; it 
differs from the result for a perfectly diffuse moon by 
a factor of 8/3. The difference between the theory 
just outlined and the smooth moon model is, of 
course, that the pow^er in each case is scattered from 
different regions of the surface as we shall see in the 
next section. 

In order to obtain a quantitative comparison of 
the results just derived with experimental data, it is 
necessary to assign a reasonable value to the power 
reflection coefficient a. If it is assumed that the 
lunar surface is composed of material with an 
effective relative permittivity €r, then the coefficient 
a can be estimated by calculating the fractional 
power scattered by a semi-infinite dielectric slab 
exposed to a normally incident plane wave: 



LVe. + lJ 



(4-19) 



In order to complete the calculation an estimate 
must be made of the relative permittivity of the 
lunar surface. 

The composition of the lunar crust has been the 
subject of considerable speculation for a great many 
years. On the basis of many studies at infrared 
wavelengths the surface appears to be composed of 
a granular rocky substratum lying beneath a very 
thin layer of dust. Thus it appears reasonable to 
conclude that for our purposes the lunar surface is 
composed of dust-like granular rock. 

Of course, little is known about the precise types 
of rocky substances to b e found in the lunar crust . On 
the other hand, many well-considered ^'guesses" 
have been made in order to provide estimates of the 
physical constants of the lunar surface. A typical 
value of the various relative permittivities of such 
substances is close to 8. Taking this to be a repre- 
sentative value for lunar surface substances in the 
solid state, an estimate can be made of the effective 
permittivit}^ of the material in a dust-like or granular 
state using formulas such as those developed by 
DeLoor [1958]. 

It can thus be shown that the eft'ective relative 
permittivity of an aggregate of closely packed spheri- 
cal granules of rock, in vacuo, whose relative per- 
mittivity €i in the solid state is 8, is equal roughly to 



4. Substituting this value into (4-19), the theoreti- 
cal radar cross section of the moon is approximately 



(7. = 0.11 ira^. 



(4-20) 



Ample data exist to check this result against experi- 
mental findings. Senior and Siegel [1959], for ex- 
ample, present a table of several measured values of 
(Ts/TTd^ over a considerable frequency range. In a 
recent summary paper on the problem of lunar radar 
reflections Evans [1961] contributes further values at 
higher frequencies. As pointed out by Evans, the 
radar cross section of the moon shows no systematic 
variation with wavelength over the interval 3 m to 
3 cm and the average value obtained experimentally 
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0-,= (0.081 ±0.005)7ra2 



(4-21) 



which agrees well with the result in (4-20). 



5. Pulse Transmission: The Quasi-Specular 
Component 

Lunar radar reflection experiments in which short 
pulses were transmitted provided somewhat more 
insight into the nature of the moon as a scatterer 
of electromagnetic radiation than could be gained 
from continuous wave transmission. In the case of 
very short transmitted pulses the temporal distribu- 
tion of echo power can be interpreted as a measure 
of effectiveness of various parts of the moon^s surface 
in scattering radiation back to the Earth. As 
already mentioned, the average power received at 
any instant can be thought of as the sum of the quasi- 
specular and the diffuse components. The quasi- 
specular component dominates that part of the 
return which originates in the vicinity of the sub- 
terrestrial point on the moon. The discussion which 
follows is devoted to the derivation of expressions 
for power received from quasi-specular scattering as 
a function of time. The scattering model is the 
same as that described in section 3. 

In the case of transmission of a rectangular pulse 
of duration r the incident magnetic field can be 
expressed in the form 

W=eiHiexp[—i{kx^+c^t)][U{ct+x^)-U{ct-CT+Xs)] 

1 r"° 

=^i^i2^ exp [i{ct+xs)ki] [1 

-exp [-i{k+h)cr]]^~^dK (5-1) 

where U denotes the Heaviside step function. It 
should be remarked that the path of integration in 
the second expression passes below the pole on the 
real axis. If JS.i{k) denotes the incident wave in 
the case of continuous wave transmission, apart 
from the time factor exp {—icot), then an alternative 
expression to (5-1) is 
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+k)cr]]Kd-ko[j~jdK (5-2) 

Thus if Hs(A:) denotes the scattered field (4-5) with- 
, out tlu^ factor exp (—io)t), then the scattered magnetic 
field in the case of pulse transmission is just 

1 r°° . 

H'=,^^ exp {ikict)[l—exp [—i{ki 

+k)cr]]H,(-k,)[j^yku (5-3) 

As in the previous analysis, the scattered electric field 
can be obtained from the relation 



E^=-foe2/^^ 



(5-4) 



These expressions will be used subsequently to cal- 
culate the expected value of received power as a 
function of time. 

Once again, the quantity of interest is the expected 
value of the Xs-component of the Poynting vector S 
at the position of the observer. In the case of 
continuous wave transmission, the rapid fluctuations 
of S were suppressed by calculating a local lime- 
averaged value. In the case of pulse reflection, the 
same time average must be interpreted as a local 
^^adiabatic" average in which the time dependent 
coefficients of the harmonic factors are regarded as 
slowly varying over several periods of the carrier 
wave. In the case of gross roughness, this is a good 
approximation except possibly in a few very narrow 
intervals in the time donniin (such as the expected 
position of the leading edge of the return). Since 
experimental data are uncertain in these regions, the 
precise time average at these points deserves no 
special consideration. Thus for pulse transmission 
the quantitv to be derived is the stochastic average 
of 



S=^ Re E^Xff*. 



(5-5) 



Just as in the case of CW transmission the present 
set of assumptions leads to a Poynting vector in the 
x^-direction at the receiver so that Sx^ will again be 
written simply as S. 

It is convenient to shift the origin of time so that 
the ^ 'average'^ leading edge of the scattered pulse 
arrives at time T=0; i.e., put ct = cT+R—2a in 
(5-3) and (5-4). Some degree of brevity is achieved 
by setting 



F{k.^ = [l-exp [-iik,+k)cr]] [^] . 



(5-6) 



Here, r is the duration of the transmitted pulse and 
ki is a transform variable. The expected value of 
S can then be written in the form 



. exp [i{cT-2a)ki-i(cT-2a)k2] • F{k,)F\h) 

• exp [i2kia cos ^j— i2fea cos 62] 

• ^{exp [i2kihi cos 6i—i2k2h2 cos ^2] 

JziJz2]ddid62d(i)id(t)2dkidk2 • (5-7) 

It is important to note that the path of integration 
in the second transform (involving ^2) passes over 
the pole on the real axis since we have taken the 
complex conjugate of the expression for H^. 

In the remainder of the discussion we shall confine 
our attention to pulses of very short duration 
(r<<<2a/c) and to comparatively small values of 
time T. With these restrictions, the same kind of 
approximations used in the CW analysis will apply 
to the present situation. For example, when T is 
very small compared with 2alc, that portion of the 
lunar hemisphere contributing to the scattered 
power is characterized by small a^^^sin di. More- 
over, the receiver bandw^idths employed in lunar 
scattering experiments are sucli that only those 
frequencies within a fc^w t (mis of megacycles of the car- 
rier frequency are of significance from the equipment 
standpoint. Although the variation of k\ and k2 in 
the transforms corresponds to all frequencies, in so 
far as the geometric approximations are concerned, 
products such as kikjo^ can be treated as being of 
tlie order of Fo-^ Hence, apart from a few obvious 
modifications, the simplification of the integrand in 
(5-7) proceeds in the same fashion as before, giving 
as a result 

«^isi =[f «] [M "'■[£ £ "" ^""' 

-icTk2] exp [-2a'(ki-k2y] 

• F{k^) F*{k2)L{ki,k2)dk,dh'l i^-S) 
where 

^(A^i,*2)= JJJ^'JJ'JJ' exp [~i (hr + \')x{^ 

•exp! -(^J'2-X')^iJ 

• exp [X^XiX2 cos (01—02)] dx\dxld4>id<j>2. (5-9) 

In this last expression, }^=^kik20^\p\^)\, vi = 2kiay 
and V2 = 2k2a. 

The integration in (5-9) is performed in nmch the 
same manner as that in (4-14). The result can be 
expressed in the form 

^^^"^^)=^[^+^'-^(^'-^^)]"' ^^-10) 

where ^=2a/\p'\0)\a\ 

The calculation of the resulting transforms in 
(5-8) is somewhat tedious and will not be reproduced 
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here. The final expression for the expectation of S, 
under current restrictions, can be conveniently 
written as 



^{s}=[- 



where 



fo 



m 



^[J(cT)-J{cT-ct)]], 



(5-11) 



J(^) = exp (-20V)[l + erf {u/2^f2a)] 

- exp i-^u) [l +erf (u^2^|2a- ^J2^(J)] . (5-12) 

Although (5-11) is a relatively uncomplicated 
relation and well suited to rapid calculations by an 
electronic digital computer, more_ insight into the 
gross functional dependence of S^'xS] on time is pro- 
vided by further simplification. We first observe 
that the product 13a is a number very much less than 
unity (except for completely negligible lunar surface 
slopes). Furthermore, on the basis of several cal- 
culations described elsewhere [Winter, 1961], it has 
been shown that the principal role of the error 
functions appearing in (5-12) is to smooth out dis- 
continuities in slope which would otherwise appear 
if a were allowed to approach 0+ in their arguments.^ 
Therefore, replacing 13a by zero and noting that 

lim ~ [l+evf{u/2^2a)] = Uiu), (5-13) 

it can be shown that an approximate relation for the 
expected value of S is 



ro 

l-exp i—pcT) 



; T<0 
; 0<T<r 



L [exp (^cr)-l] exp (-^cT); T<T«2a/c. 

(5-14) 

Thus for values of T greater than the pulse duration 
the decay of echo power is approximately exponential, 
with a ''time constant'' equal to 1/^c, 

6. Pulse Transmission: The Diffuse 
Component 

The quasi-specular component of the lunar echo 
is due to scattering from that portion of the lunar 
terrain which appears to be locally smooth and 
whose irregularities are large compared with a wave- 
length. In the case of short pulse transmission this 
component will decrease in importance very rapidly 
as the '^illuminated zone'^ moves back toward the 
limb from the subterrestrial point. At the same 
time scattering from smaller, rougher surface features 



^ This is equivalent to letting a->0-\- in the factor exp [— 2a2(A:i— /;2)2] appearing 
in (5-8). 



(such as rocks and jagged facets) will increase in 
relative importance since these facets will tend to 
scatter radiation in a more isotropic fashion. Since 
the detailed structure of such features is not known, 
a semiempirical approach must be employed. 

We begin by assuming that echo power associated 
with the diffuse component is due to scattering from 
a great number of small scatterers which are ran- 
domly distributed over the surface of the moon. 
Suppose that the jth scatterer at (dj,<j)j) is situated 
within the small area element AAj. The average 
total power incident upon AAj is approximately 

P,„e = [f ^l]/(/?,) COS dAAj, (6-1) 

where Rj is the range from the transmitter-receiver ^ 
to the scatterer and the factor /(/?;) is such that if 
g[cT—{Ej + a—R)] is the instantaneous distribution 
of transmitted power in range, then 



f{R,)=g[cT-2{Rj+a-B)]. 



(6-2) 



Sj=P.-^- 



'^'"'liwRU 



<^kiS),4>j) sec Bj. 



(6-4) 



Thus for a rectangular pulse of duration r arriving j 
at the subterrestrial point on the moon at time 
T=0, the factor g takes the form 

g[cT-{Ej+a-R)]=U[cT-{Rj+a-B)] 

-U[cT-CT-iRj+a-R)]. (6-3) ^ 

Let 'ak(dj,(t>j) denote the radar cross section per unit 
area of the ^th scatterer. The power flux at the 
receiver due to scattering from this facet is then 



Since it has been assumed that the scatterers are ^ 
distributed over the lunar surface in a random 
fashion, the relative phases of waves reflected from 
them will be randomly distributed and uncorrelated 
and the power flux at the receiver at any instant of 
time is simply the sum of the individual contribu- 
tions. Combining (6-1), (6-2), and (6-3) with (6-4) 
the received power flux can be written as 

^{S}^h{Jt> {UlcT-2{Rj+a-R)] 
-U[cT-CT-2{Rj+a-R)]] 

'\j^^M9j,cl>j)AAj, (6-5) 

where Njs the number of scatterers. 

Since ajc{dj,4>j) is not Imown, it is convenient to 
treat it as a random variable with a probability 
density function i?pA;(^jf,0j)]. The expected value of 
S may then be defined as 

^ {S} =S S(a,)p[a,(ej,h)]Aa,. (6-6) 
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Equation (6-6) can be coinbinecl with ((j-d) to give 
an explicit expression for rf{S}. If we allow the 
number of scatterers to approach infinity and let all 
A(7/, to go to zero (thus forcing M to approach infinitv) 
the sunnnations become integrations over 6, cj), and a. 
The range Rj is then given by R~a cos 6. The factor 
ll4tTrR] may be replaced by l/47r7?^ and subsequently 
be removed from the integral. With these modifica- 
tions, the expectation of S takes the form 



■^'^'=[i-o[.-i]rr,r 



\U{cT-2a+2a cos &)-U{cT~cr-2a+2a cos ^)1 

. a (^,</)) p [^ ((9,0) ] sin ddddctxia. (6-7) 

The integral over a is nothing more than the mean 
a^o(^,</>) of the random variable o{d,(t>). Thus, in order 
to complete the calculation some additional state- 
ments must be made with regard to o^o(^,0). 

As was stated earlier, in the absence of any detailed 
information concerning the structure of the less direc- 
tive scattering facets, it is impossible to assign a pre- 
cise functional form to ao{d,<t>)' On the other hand, 
since the facets are presumably distributed over the 
surface in a random fashion, there is no reason to 
retain an azimuthal dependence explicitly. This de- 
pendence could, in fact, be determined (in theory) 
only if the configuration of the facets were specified 
together with their orientation with respect to the 
direction of the incident field vectors. In suppress- 
ing the azinnitlial depend(Mice, tlie parameter a„(^) 
can be interpreted as the vahu^ of ao(^,(/)) averaged 
over the angle (^. 

The experimentally determined behavior of_tlie 
diffuse component (e.g., see fig. 1) suggests that ao{d) 
may be represented in the form 



C7o(<9) = o-o(cos6>)", 



(6-8) 



where n is of the order of magnitude of unity or juore. 
If the surface facets reradiated energy isotropically, 
and shadowing were unhnportant, then n would be 
identically equal to unity. On the other hand, we 
shall see in section 7 that PettengilPs experiments 
suggest the net effect of backscattering from these 
facets is more adequately described b}^ Lambert ^s 
Law which implies that n is closer to two. It should 
not be concluded from this, however, that the fea- 
tures of the lunar surface producing the diffuse echo 
component scatter radiation in the same sense that a 
'^gray body" is diffuse at infrared wavelengths. In 
fact, if further experimentation verifies this law, and 
if the present model is correct, the diffuse component 
may be interpreted as scattering from certain pre- 
ferred types of configurations, such as spikes or other 
sliarply pointed facets. Although the following is 
not particularly appropriate to the physical situation 
at hand. Anient [1956] has shown, by way of example, 
that Lambert's Law is suggested in backscattering 
from the coplanar ends of parallel but randomly situ- 
ated semi-infinite thin wires. 



At any rate, when (6-8) is substituted into (6-7), 
the scattered power flux associated with the diffuse 
component can be written as 



.f{s}==[|^?][^T ^^^ 



_2RJ 71+1 



'< 



[[■ 






; T<o 



0<T<r 



r<T< 



2a 



c 

(6-9) 



over the time interval of interest here. The total 
lunar echo power at any instant of time 7" is the sum 
of (6-9) and the corresponding expression for (5''{S} 
associated with quasi-specular reflection. Tlie value 
of Co is, of course, determined from (\\j)eT'im(Mital 
data. 

7. Comparison of Total Theoretical Signal 
With Experiment 

Altbough experimental attempts to obtain lunar 
echoes using short pulse transmission have been in 
progress throughout tlie past decade, detailed results 
of such studies have appeared in the literature only 
within three or four years of this writing. Inade- 
(juacies in existing ecpiipnuMit, coupled with large 
transmission losses, still make it difficult to attach a 
great amount of accuracy to any collection of lunar 
radar observations. Therefore, it is not unadvisable 
to regard the experimental results available at the 
present time with some reservation. On the other 
hand, a number of experiments, witfi various carrier 
frequencies and pulse lengths, are sufficiently relia- 
ble so that a satisfactory degree of accuracy may be 
attributed to the measurements. Amongst these 
data are those reported by Trexler [1958], Hughes 
[1960], and Pettengih [1960, 1961]. Theoretical re- 
sults derived in the present investigation are com- 
pared with the aforementioned experimental results 
and others in the discussion below. 

Among the most accurate lunar radar data avail- 
able are those of Pettengill [1961]. The equipment 
employed in this experiment was the Millstone 
Radar of Massachusetts Institute of Technology's 
Lincoln Laboratory. Using a transmitted carrier 
frequency" of 440 Mc/s and a pulse length of 100 
jLisec, Pettengill was able to obtain an average 
temporal distribution of echo power throughout a 
time interval corresponding to the total radar 
^'depth'^ of the moon (11.6 msec). Figure 1 pre- 
sents the results as reported by Pettengill. The 
precise behavior of the initial rise of the signal from 
the vicinity of T=0 is not well known so that only 
the period of signal decay will be compared with 
theory. 

The task of fitting the theoretical results to the 
experimental curve involves the solution of a three 



(525829 — f)2- 
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parameter problem since the total signal is given 
by the sum of (5-14) and (6-9) which together 
involve the constants ^, ao, and n. Instead of 
assigning values to the constant /5 we shall refer to 
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Figure 3. Comparison between theory and experiinent ove?- 
an entire echo. 
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Figure 4. Comparison of theoretical results with the initial 
portion of the echo decay for a very short pulse. 



the combination 1/2 Vai^- From the definition of 13 
and the remarks in section 4 with regard to the 
relation between |p''(0)| and the mean square lunar 
surface slope, it can be seen that 



m<m""'o^^HHi^mi 



(7-1) 

In the special case of a Gaussian correlation 
function, the various expressions in (7-1) are equal 
simply to af8 (for 5<<a). 

On the basis of several calculations it was con- 
cluded that the best overall agreement between the 
theoretical curve and the experimental one is 
achieved by using parameter values of the order of 

1/15 <l/2VttiS< 1/12 



1/4 < 0-0 < 1/3 

3/2<7i<2. 



(7-2) 



The graph corresponding to the set (1/12, 0.25, 2) is 
shown together with the experimental decay rate in 
figure 3. 

Of course, comparison of the theory with a single 
experiment is inconclusive. The general usefulness 
of the present formulation will be demonstrated if it 
can be shown that comparisons with other experi- 
mental results lead to values for l/2-y/aP^ o-q, and n 
which are close to those indicated in (7-2). 

Recent work reported by Hughes [I960] provides 
us with results of an experiment performed with 
equipment parameters which are considerably dif- 
ferent from those of the Millstone Radar. Using a 
5 jLtsec pulse and a carrier frequenc}^ of 3000 Mc/s, 
Hughes observed a lunar echo decay which seemed 
to be fairly well represented by the function 



S=Soexp{-10df); 3°<^^<14° 



(7-3) 



where ^/=arc cos (l — cT/2a). This expression is 
plotted in figure 4. The theoretical result cor- 
responding to the set of parameter values {l|2^jafi, 
o'o, ^) = (1/15, 0.25, 2) appears to give a fairly good 
overall fit to the data. 

Inspection of figures 3 and 4 reveals a s^^stematic 
difference in slope between the theoretical and 
experimental results even though the parameter 
values fall within the ranges indicated in (7-2). 
Similar disparity in slope was encountered by 
Hagfors [1961] who assumed a Gaussian correlation 
function to derive an ^^apparent delay spectrum, '' 
which can be interpreted as a measure of the relative 
effectiveness of various portions of the lunar surface 
in returning a radar signal. Hagfors remarks that 
the discrepancy in slope might be removed by a 
different choice for the correlation function p. 
Evans [1961] concurs with this view. Another 
possible explanation is suggested by an inspection 
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of luiiai' |)l)Olo<;Taphs or charts of lunar topography 
which reveals that the subterrestrial ponit hes within 
a comparatively smooth flat region which gradually 
bl(Mi(ls into somewhat rougher terrain around its 
borch^rs. Thus it uui}^ be that the effective rms 
value of local lunar surface slope actually mcreases 
slightly with increasing d in this general region. In 
this event the predicted rate of decays of echo power 
with time would slowly decrease with increasing T 
as see] US to be observed in practice. 

'J'heories based on the statistical model of the 
moon's surface have been criticized in the past for 
their inability to account for the so-called modulation 
loss phenomenon. The term ^'modulation loss" re- 
fers to the apparent decrease in the peak power of 
the lunar echo as the transmitted pulse duration is 
diminished. The present formulation, on the other 
hand, |)rovides an explicit expression for the modu- 
lation loss for comparison with observed data. Al- 
though the appropriate experimental results are 
sparse and somewhat uncertain, the most reliabh* 
are probably those communicated to Senior and 
Siegel [1960] by Youmans. Acconhng to these data 
the peak |)ower returned from a 200 /xsec pulse is 
about 1 (lb ])elow the value for continuous wave 
transmission. The corresj^onding value for a 30 
jusec pulse is 8.5 clb, for a 10 ^sec ])ulse it is 17 db, 
and for a 2 /xsec pulse the loss appears to be 22 ±2 
db. An approximate theoretical expression for the 
modulation loss as a function of pulse length is 
obtained from a combination of (4-17) and (5-14): 



Loss=l — exp (—I3ct) 



(7-4) 



Figure 5 presents a plot of the modulalion loss 
predicted from (7-4) with l/2VaiS=l/12 togethcM- with 
the experimental results described above. The 
agreement is seen to be quite satisfactory. 

Another check of the foregoing results was per- 
formed using the data obtained by Trexler [195<S| 
who employed a pulse length of 12 jusec and a carrier 
frequency of 19(S AIc/s. Trexler presents a plot of 
the fraction of total echo energy received at time T. 
With (5-14) it is a simple matter to obtain an expres- 
sion for this fraction as a function of time. De- 
noting this quantity" by Q, for T> t, we have 

Q=l-{^cr)-' [exp (^cr)-l] exp (-^cT). (7-5) 

Figure 6 displays curves of log^ (1 — Q) as a function of 
T for T=12 jLtsec. The circled points were calcu- 
lated from Trexler's plot of experimental results. 
It would appear that the best fit to the radar data is 
obtained by using the ratio 1/2Vg^/5= 1/17.5, a value 
slightly lower than indicated in (7-2)^^ The the- 
oretical curve corresponding to l/2^lal3=l/15 is 
shown for comparison. 

A check of tlie paranu^ters a^ niu] n is not readily 
available since oidy one or two lunar radar equip- 
ments have sufficient sensitivity to detect the low 
level deca}' associated with the tail of the echo. 
However, Pettengill [1960] has also reported the 
results of experiments with the Millstone Radar in 



u 


1 


1 1 — 


. 


1 

n 


-•s 


// 




1 


1 




2v 


/^ 


12 




r 








10 


' 


THEORY 

O EXPERIMENT 






15 


\-l 









Figure 5. 



50 100 150 200 250 

PULSE LENGTH r./usec 

Comparison between theoretical and experimental 
modulation loss. 






o 


1 




1 
1 1 






o ^\ 


^^-^ 


^ 2y^ '5 


-1 




o 


\^ 


^^^\ 


-2 


- 


1 


1 
17.5 


^^^ ^s. O 


2^0^ 


-3 


— 


THEORY 








o 


EXPERIMENT 












T= 12 /j.sec 


-4 




-^ 


L 


i _ . 



50 



100 
TIME T ,fxsec 



200 



Figure 6. Comparison between theoretical and observed 
fraction of echo energy returned at time T. 



which 500 jLisec pulses were used. That portion of 
the experimental results corresponding to the diffuse 
component is illustrated in figure 7 togetherjwith 
the theoretical curve corresponding to (l/2-\ial3, o-q, 
n) = {l/l2, 0.3, 2). Agahi the parameters fall within 
the ranges indicated by (7-2) and overall agreement 
is quite good. Comparison with portions of the 
experimental curve corresponding to times less than 
1 msec is not possible because the sampling interval 
employed in the experiment (500 jusec) is so great as 
to obscure the details of the rise and decay of the 
echo in this region. 
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Figure 7. Comparison of theoretical and experimental results 
for the diffuse component. 
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